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Abstract 

ra-DBI gravity is a gravitational theory which yields near de Sitter inflation sponta- 
neously at the cost of breaking Lorentz invariance by a preferred choice of foliation. We 
show that this breakdown endows ra-DBI gravity with one extra physical gravitational 
degree of freedom: a scalar graviton. Its existence is established by Dirac's theory of 
constrained systems. Firstly, studying scalar perturbations around Minkowski space- 
time, we show that there exists one scalar degree of freedom and identify it in terms of 
the metric perturbations. Then, a general analysis is made in the canonical formalism, 
using ADM variables. It is useful to introduce an auxiliary scalar field, which allows 
recasting ra-DBI gravity in an Einstein-Hilbert form but in a Jordan frame. Identifying 
the constraints and their classes we confirm the existence of an extra degree of freedom 
in the full theory, besides the two usual tensorial modes of the graviton. We then argue 
that, unlike the case of (the original proposal for) Hofava-Lifschitz gravity, there is no 
evidence that the extra degree of freedom originates pathologies, such as vanishing 
lapse, instabilities and strong self-coupling at low energy scales. 
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1 Introduction 

General diffeomorphism invariance is a central property of General Relativity (GR). In par- 
ticular, it means that all foliations of space-time (say by space-like hypersurfaces) are equiv- 
alent. For many pratical purposes a choice of foliation must be made, which can be achieved, 
for instance, by an ADM decomposition of space-time [HI H] . The ten metric degrees of free- 
dom are then decomposed into the lapse, shift and 3-metric, of which only the latter contains 
physical information. The lapse and the shift are pure gauge, as it is manifest in the fact that 
their equations of motion are actually constraints. Such constraints need only be imposed 
in one leaf of the foliation and are guaranteed to be preserved in the evolution. The number 
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of physical degrees of freedom of the theory can be formally computed by Dirac's theory of 
constrained systems [5]. For GR there are 20 phase space variables and 8 first class con- 
straints, leading to (20 — 2 x 8)/2 degrees of freedom. These are the two polarisations of the 
graviton, which are clearly seen by studying plane waves in linearised theory, corresponding 
to trace-free, transverse, tensor modes. 

Any theory of gravity for which general diffeomorphism invariance is broken will have, 
generically, a different number of physical degrees of freedom. An example is when the full 
diffeomorphism invariance is broken to foliation preserving diffeormorphisms (FPD). This 
type of symmetry has been extensively discussed in recent years in the context of Hofava- 
Lifschitz (HL) gravity [6]. Under this smaller symmetry, the shift can still be gauged away, 
but the lapse (and the conformal mode, cf. Section 13.1.21) acquires physical significance. As 
a manifestation, the equation of motion of the lapse, the Hamiltonian constraint "H, is not 
preserved in the evolution: the evolution of H, is actually a dynamical equation for the lapse. 
Naively, this seems to introduce one new degree of freedom. In the non-projectable class of 
HL theories, however, a formal counting, based on the fact that the Hamiltonian constraint 
becomes second class and the lapse is determined (removing two phase space variables) , gives 
(2x (6 + 3)— 2x (3 + 3) — l)/2 = 5/2 degrees of freedom [10]. Getting one half degree of 
freedom might not necessarily be problematic in itselfQ There are, however, at least three 
problematic features related to the existence of this new half degree of freedom in HL gravity. 

The first problem concerns the absence of dynamics and is associated to the fact that the 
extra degree of freedom is halved. For generic asymptotically flat space-times, the lapse is 
forced to be zero at spatial infinity [10J. Actually, for particular values of the couplings, it was 
shown (and suggested this could be the case for generic values of the couplings) that the lapse 
must vanish everywhere. This indicates that there is no dynamics in HL gravity. The second 
problem is the short distance instability the scalar mode might trigger (if the dynamics was 
not frozen). Looking at perturbations around generic backgrounds, it was found in [IT] that 
the high frequency modes of the extra degree of freedom develop an imaginary part and the 
perturbations can grow very swiftly in tinted The third problem is the self-coupling of the 
scalar mode which remains strong to very low energy scales [Til IT2] . This implies that the 
extra scalar mode never decouples and thus HL gravity does not flow to GR in the IR as it 
was hoped. 

Recently, another model of gravity for which the full diffeomorphism invariance is broken 
to FPD was proposed by some of us [HE]: n-DBI gravity. This model was designed to 
reduce to the Dirac-Born-Infeld (DBI) scalar field theory for conformally flat geometries; 
it was motivated by the observation that for such geometries, in a cosmological setup, two 
epochs of acceleration for the universe are naturally obtained [TJ . Subsequently, we observed 
in [2] that a large class of exact solutions of GR (including the most standard spherical black 

*A well-known example is a chiral boson. 

2 The first two issues are apparently contradictory; if the lapse must collapse everywhere, there would 
be no way to develop an exponentially growing mode involving the lapse. We will, however, discuss the 
instabilities of the model which evades the first issue. 



3 



hole solutions, in the most standard foliations) are solutions of n-DBI gravity. 

The purpose of the present paper is twofold. Firstly, to understand what are the extra 
degrees of freedom present in n-DBI gravity due to the breakdown of full diffeomorphism 
invariance. Secondly, to investigate whether n-DBI gravity is afflicted by problems similar 
to those described above in the context of HL gravity. 

Concerning the first part of our investigation, we shall show that n-DBI gravity pos- 
sesses one extra degree of freedom, as compared to GR, that this degree of freedom is a 
scalar graviton and that its existence is a direct consequence of the breakdown of the full 
diffeomorphisms group to FPD. This shall be done in Sections 3 and 4, first by considering 
perturbations around Minkowski spacetime and then by considering the fully non-linear n- 
DBI theory. Concerning the second part of our investigation, we shall present arguments, in 
Section 5, that none of the aforementioned problems that afflict HL gravity are present in 
n-DBI. This healthy behaviour is intimately connected with the nonlinear lapse dependence 
in n-DBI, as compared to HL gravity. Indeed, n-DBI gravity naturally contains lapse de- 
pendent terms in the action, analogous to those suggested in [7] for the consistent extension 
of HL gravity. 

The analysis in Sections 4 and 5 can be made simpler by introducing alternative forms for 
n-DBI gravity. One alternative form linearizes the action, i.e. eliminates the awkward square 
root dependence, by introducing an auxiliary field. Another alternative form makes the 
action covariant, by introducing a Stiickelberg field. These alternative forms are introduced 
in Section 2. We shall conclude, in Section 6, by discussing more on the stability and 
commenting on future research directions for understanding this model. Some technicalities 
are presented in three Appendices. 



2 Different forms of n-DBI gravity 

The action for n-DBI gravity without matter is [H [2] 

where Gat is Newton's constant, q, A are the two dimensionless parameters of the theory and 
R is the four dimensional Ricci scalar. To completely define the theory a foliation structure 
must be chosen. Let n be a unit time-like vector field, everywhere orthogonal to the leaves 
of such foliation; let h^ u = g^ u + n M n„ be the first fundamental form and K^ u = \£ xl h^ v the 
second fundamental form. Then /C in (12.11) is defined as 

/C = — -=£ a (VhK) , K ee K^hT . (2.2) 
yh 

Performing the Arnowitt-Deser-Misner (ADM) decomposition [3] 

ds 2 = -N 2 dt 2 + hij (dx i + N*dt) (dx j + N j dt) , (2.3) 
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where N and Ni, respectively, are the lapse and shift functions, and hij with i,j 
the spatial metric, the action (12. ip becomes 

3A 



S 



dtd 3 xVhN 



1 + ^(R + KuK* 



K 2 - 2N- X AN) - q 



1,2,3 is 



(2.4) 



where R is the Ricci scalar of h^. In particular, observe the lapse term N~ l AN. As suggested 
in [TJ, adding analogous terms can provide a consistent extension of HL gravity. Indeed, as 
we shall see below, this term plays a crucial role in evading the pathologies that afflicted HL 
gravity. 



2.1 A linearised form with auxiliary field 

As it is known in various contexts, square root type actions may be linearised by introducing 
auxiliary fields. A well known example is the classical equivalence between the Polyakov and 
Nambu-Goto actions in string theory, via the introduction of an auxiliary metric (the world 
sheet metric). A similar reformulation for the Eddington inspired Born-Infeld gravity was 
provided recently [8J, for which the auxiliary variable is the "apparent" metric. In our case, 
with the introduction of an auxiliary scalar field e, the action (12.11) becomes^ 

S e = ~\c~~r i / <*V=0e \ W R - 2G^A c (e) + /c] , A c (e) =®±-(*l-l-L) . 
16ii& N J L J G Z N V e e 1 ) 

(2.5) 

For constant field e, this is the Einstein-Hilbert action with a cosmological constant (in 
exactly the form found in [2] in terms of the integration constant C therein) with the Gibbons- 
Hawking- York boundary term [Hill]. For generic e, the theory resembles General Relativity 
in a Jordan frame (but indeed it is quite different). This form will be used in Section H] for 
the canonical formalism analysis of n-DBI gravity. 

It is worth noting that we can rewrite the action (12.51) in the Einstein frame by performing 
a Weyl transformation: 

gp,v -»> e^g^ . (2.6) 
Redefining the auxiliary field e = exp (2%), the action becomes, up to boundary terms, 

Einstein = "Tr^ / <*S/=0 - 6 + h a P) d aX dp X + + V( X )] , (2.7) 

lb7rGAr J L J 

with the potential 

6A 

^(x) = 7^exp(-4 X )[(exp(x)-exp(-x)) 2 + 2(l-g)] . (2.8) 

There is a caveat: it is misleading to regard the Einstein frame theory as a scalar-tensor 
theory (the scalar being x). Despite the appearance of the kinetic term, the scalar field x is 
still an auxiliary field and does not give rise to an independent degree of freedom. As we will 
see, the extra scalar mode is furnished in the metric and the scalar field x is only related to 
it through the equations of motion. 



3 We thank Soo-Jong Rey and Takao Suyama for suggesting this formulation. 
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2.2 A covariant form with Stiickelberg field 



n-DBI gravity breaks Lorentz invariance due to the coupling of gravity to the unit time-like 
vector field n, which defines a preferred space-time foliation. Full general covariance can 



be restored by introducing a Stiickelberg field 
time-like and non- vanishing: 



0„ 



n, 



(x^ 



-1 



such that its gradient is everywhere 



x = <r a, 



(2.9) 



By definition the theory is invariant under — » /(</>), where /(</>) is an arbitrary function of 0. 
Note that the original non-covariant form is recovered when = t for which n = (— N, 0, 0, 0), 
and the symmetry reduces to the time reparametrization t — > fit). A similar treatment in 
the case of HL gravity has been performed in [llj. The extrinsic curvature becomes 



K 



[IV 



[n a D a h^ u + h m D v n a + h ua D^n c 



K = D n n a 



It then follows that 

K = -2D a (n a D^) . 
Thus, the n-DBI action may be rewritten in the covariant form 



3A 



d x 




N /(4) 



U - 2D C 



-,D 



dp<t> 



(2.10) 
(2.11) 

(2.12) 



This form will be used in the discussions of Section |5j In contrast to the Einstein frame 
theory, this covariant theory can be thought of as a scalar-tensor theory. Indeed, in this 
form, the scalar field does yield an independent degree of freedom. Put differently, the 
scalar mode in the metric is entirely transferred to the Stiickelberg field <fi. 
Clearly, this can be linearised again by introducing the auxiliary field e, 



St 



16ttG 



jV 



(fx^ge l W R - 2G N A c {e) - 2D C 



--D, 



(2.13) 



In these formulations, general covariance gets spontaneously broken by the scalar field 
acquiring the vev (0) = t. It should be noted that this has a certain bearing on the ghost 
condensation of |13|. 



3 The scalar graviton 

In this section we shall show that there are three physical degrees of freedom in n-DBI 
gravity: the two usual tensorial modes of general relativity plus one "scalar graviton" . This 
is in contrast with HL gravity where there are two and a half physical degrees of freedom 

nm. 
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3.1 Scalar perturbations around flat space-time 



We shall first study the scalar mode in flat space-time. In particular, we shall clarify how the 
scalar mode becomes physical as a direct consequence of the breakdown of full diffeomorphism 
invariance to FPD. 



3.1.1 The existence of a scalar graviton 

The metric has ten components; the lapse N, the shift iV l , and the spatial metric h^. These 
can be decomposed into 4 scalars (n, B,ip, E), 2 transverse vectors (Ai,Ai), 1 transverse 
traceless tensor hij as follows (see e.g. [H]): 

N=l + n, Ni = ViB + Ai , 

=6ij - 2 (6^ - t/j - 2^^£ + (Viij + VjA^j + hij , (3.1) 

where the transversality and traceless conditions are imposed, 

V*A = V l Ai = V%j = ~h l i = . (3.2) 

We now expand the n-DBI gravity action (12 .4p to quadratic order around flat space- 
time, setting q = 1. Thanks to conditions ( 13.21) . scalar perturbations and vector-tensor 
perturbations decouple from each other. The quadratic Lagrangian density for the scalar 
fields yields, up to total divergences, 

47rG;v£ scalar = 2ip 2 + Aip(E + Ab) + (2ip - An) Aip + — ^ (2A^ - An) 2 . (3.3) 

V / 6A 

This essentially corresponds to the A = 1 case of [TJ, and the scalar graviton in this model is 
qualitatively different in its character from the one discussed there. 

In order to establish the existence of the scalar degree of freedom, we shall apply Dirac's 
theory of constrained systems. The Hamiltonian density of this effective scalar theory read^l 

Scalar = + l*** ~ " ^ ~ ^ ^ ~ if ^ ~ ^ ' (3 " 4) 

where pe and p^ are the conjugate momenta of E and ip, respectively. Since the Lagrangian 
density f 13 . 3 jl does not contain the time derivative of n or B, their conjugate momenta become 
the following primary constraints: 

$ 1 =p n = 0, <^ 2 =p B = 0. (3.5) 

Thus the dynamics of this system is governed by the Hamiltonian density 

Scalar = Scalar + + ^2 , (3.6) 



'For convenience, we have rescaled the scalar fields by the factor of (AttGn) 1 ^ 2 ■ 
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where Ai and A2 are Lagrange multipliers. The consistency requires the time flows of the 
primary constraints to vanish (herein « means 'weakly equal' as standard in Dirac's theory): 



$ 4 = $ x = 4A^ + (4A 2 ^ - 2A 2 n) « , $ 5 = $ 2 = -Ap £ « . 



(3.7) 



The former corresponds to the Hamiltonian constraint and the latter to the momentum 
constraint, as will become clear in the next section. The time flows of these secondary 
constraints do not yield new constraints. 

Among these four constraints, $2 and $5 commute with all the constraints and are 
therefore first class, whereas $1 and $4 are second class. Hence the physical degrees of 
freedom of this system are counted as (2 x 4 — 2 x 2 — 2)/2 = 1. In contrast, in GR where 
A —7- 00, all the constraints are first class and the counting becomes (2x4 — 2x2 — 2 x 
2)/2 = 0. In other words, GR has two pairs of first class constraints associated with the 
full diffeomorphism, whereas ra-DBI has only one pair reflecting that only FPD is preserved. 
Thus, in n-DBI gravity there remains one physical scalar degree of freedom that cannot be 
gauged away. We will elaborate on this point in a moment. 

Note that it is the presence of the lapse term (An) 2 ~ (N^ 1 AN) 2 in the Hamiltonian 
that leaves one physical scalar degree of freedom, as opposed to 1/2 in HL gravity^ This is 
similar to the consistent extension of HL gravity [7] . 

3.1.2 The identity of scalar graviton 

Having established the existence of the scalar graviton, we wish to identify it more explicitly 
in terms of the metric perturbations. For this purpose, we first note that the scalar field 
theory fl3.3[) has the FPD gauge symmetry: 



where T and L are related to infinitesimal coordinate transformations by £° = T and £ l = 
V J L. In n-DBI gravity, T is a function of time only, whereas L is a function of both space 

In non-projectable HL gravity, the lapse N is a Lagrange multiplier as in GR, but the Hamiltonian 
constraint is second class. Moreover, the time flow of the Hamiltonian constraint yields an additional 
constraint that depends on the lapse N. Thus, there are 3 second class constraints; the conjugate momentum 
Pn of the lapse, the Hamiltonian constraint H, and its time flow %. Together with 6 first class constraints, 
the number of degrees of freedom is counted as (2 x 10 — 2 x 6 — 3)/2 = 2 + 1/2. Note, however, that linear 
perturbations about flat space-time yield a misleading result. There appear four second and two first class 
constraints, implying incorrectly that the number of scalar degrees of freedom is zero. In projectable HL 
gravity, there is an additional primary constraint d^N = on top of p^ = 0. In this case, the time flow of 
Pn = does not yield the Hamiltonian constraint. Instead, it determines the Lagrange multiplier of diN . 
Hence pn and diN are the only second class constraints, and the total physical degrees of freedom is 2 + 1. 



B B + L-T 
E -> E - AL , 
n — > n + T , 



(3.8) 
(3.9) 
(3.10) 
(3.11) 
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and time. In GR, T also becomes a function of both space and time. Note that the gauge 
invariant quantities are 

{tp, E + AB + An} for both GR and ra-DBI , (3.12) 

{E + AB, din} for n-DBl only . (3.13) 

As it is clear from the above counting of degrees of freedom, the extra scalar graviton exists 
in ra-DBI gravity, because the constraints <£>i and $4 are second class, whereas in GR they 
are first class and generate the gauge transformations 



d 6 y( 1 (y)$ 1 (y),n(x)j = (i(x) , jy d*yUv)*MM x )] = -^(x) . (3.14) 

Comparing this with the FPD (13.8p - (13.11j) . we find that 

Ci = T, C 2 = T, (3.15) 

where we used = + E + AB). In n-DBI gravity, as we stressed above, T is a 
function of time only. Accordingly, the constraints $1 and $4 become second class and are 
not considered as generators of gauge transformations. This implies that the scalar graviton 
should involve the lapse and/or the shift whose gauge transformations are generated by T. 
This nicely fits the expectation that the scalar graviton must have something to do with the 
foliation structure which is specified by the lapse and the shift. Indeed, the equations of 
motion, 

i) = , (3.16) 
Atp = , (3.17) 
E + AB + An + = , (3.18) 



a ^ = ~6tU ~ ^' ' (3,19) 



have in the E = gauge the general solutiorg 



B(t, x) =B (x) + B 1 {x)t , (3.20) 
n(t,x) =-B x {x) -VoO) > ( 3 - 21 ) 
i/}(t,x) =Mx) , (3.22) 



with ipo( x ) related to Bi(x) by 



A more common gauge is to set the shift B = 0. We can go from the E = to the B = gauge by 
choosing the gauge parameter L(t,x) = —Bo(x)t — ^Bi(x)t 2 . This yields the conformal mode E(t,x) = 
ABo(x)t+ ^ABi(x)t 2 . Note that in either gauge the lapse n alone only accounts for a half degree of freedom 
of the scalar graviton. 
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Note that we have imposed the boundary condition that all the fields must fall off at spatial 
infinity. In other words, the functions f(x)'s appearing in the most general solution and 
obeying the Laplace equation Af(x) = are unphysical and set to zero. 

In the Hamiltonian system, a degree of freedom is the freedom to choose a pair of initial 
data for the time evolution in the phase space. We have found exactly one degree of freedom, 
i.e., the initial data specified by a pair of arbitrary functions of space, (B (x), Bi(x)). In 
GR, these could have been gauged away by choosing the gauge parameter T(t,x) = B(t,x) 
(and the Hamiltonian constraint would have enforced ipo{x) — 0). Put differently, in n-DBI 
gravity the scalar mode is the broken gauge degree of freedom T(t, x) which obeys, by taking 
the time derivative of the the Hamiltonian constraint (I3.19p . 

A 2 T(t,x) = 0, (3.24) 

where we have used n(t,x) = —T(t,x) — if>o(x) and if)(t,x) = ipo{x). A few remarks are in 
order: firstly, this is a key equation, despite its extremely simple appearance, and notably 
originates from the nonlinear lapse term (An) 2 ~ (N' 1 AN) 2 in the Hamiltonian; secondly, 
as we will see in Section 15.21 the Stiickelberg field satisfies exactly the same equation and 
thus can be identified with the broken gauge degree of freedom T(t,x). 

Observe that the frequency of the scalar mode is u = 0, and thus this is more a zero mode 
than a propagating particle mode. Nonetheless, this is the physical degree of freedom of our 
Hamiltonian system. We will postpone the interpretation of this result for the moment and 
discuss it in Section 



4 The canonical formalism — the full theory 

We shall now consider the Hamiltonian formulation of the full n-DBI gravity theory and 
confirm the existence of one extra degree of freedom. The n-DBI Lagrangian, L n DBi, is given 
by equation (12.41) . The square root makes the analysis cumbersome, and we find it more 
convenient to work in the linearised form with the auxiliary field e discussed in Section 12.11 

^nDBi = ~\ J d 3 xy^e (R + K t] K ij - K 2 - 2N~ 1 AN - 2G N A c (e)) , (4.1) 

where k = 16ttGn- In passing to the Hamiltonian formalism, we set the notation for the 
canonical conjugate momenta as 

-^nDBI ( {hiji hij 

),(N,N),(N i ,N t ),(e,e)) ^ ^ DBI ((^,^), (iV, PiV ), (iV,,^), ( e ,p e )) . 

(4.2) 

However, the time derivatives of N and iVj are absent as in GR, and so is the time derivative 
of e. Thus we have the primary constraints, 

$!=^ = 0, $* 2 =^ = 0, $3= Pe = 0. (4.3) 
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Denoting the Lagrangian density by £, the Hamiltonian density is given by 

51 — y 'hj ^nDBI 

^ 1 



VhN~ 



2 _ ,,\ \fh~N\ k 2 



^iP 13 + 



^ [p^Pij ~^p 2 )+e(R- 2G N A c (e)) 



-Vh(eAN) 

K 



where 



P 



t 3 _ °-^nDBI 
Skin 



(K ij - h ij K)e 



The time flow of the constraints are generated by the extended Hamiltonian density 



H 



e(l) _ ,T,e(0) 



DBI ~~ DBI 



^DBI + + ^2 + A 3 *3 , 



(4.4) 



(4.5) 



(4.6) 



where Ai, A 2 i and A3 are the Lagrange multipliers. Thus the primary constraints evolve in 
time as 



a=l,3 



+ E / rfM$a(y),$i(x)}A a + /^ 2 (?/)^iW}A 



^ (p% -^ 2 ) < (7?-2G'. v A ( .(0)-2A ( 



= $4 UP J 



(4.7) 



a=l,3 



2Vy j = ¥ 5 (x) 



(4.8) 



$,(x) 



rfM^D ( B ) i(2/)^3(^)}+ J2 / A{^(2/),^3(^)}A a + / dM^ 2 (2/)^3^)}A 2i 

a=l,3^ 



+ AN = %(x) . 



~ [p ij Pij -^ P 2 )-e(R- 2G N A c (e)) 



12A 

G N V e 



(4.9) 



Therefore, in addition to the primary constraints (14. 3p . we have the secondary constraints 

$4=p N &0, $ 4 5 =^^0, $ 6 =p e ^0. (4.10) 

Finally, the time flows of the secondary constraints do not yield any further constraints, as 
shown in Appendix 1X1 
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4.1 Constraints and counting degrees of freedom 

It is easy to understand the physical meaning of the above constraints. Firstly, we can solve 
$6 = for e and obtain 

e = ±(l + ^nY 1/2 , (4.11) 



6A 

where we defined 1Z = + /C. Choosing the positive sign and plugging it into $ 4 , we find 



^ K 
*4 = -2 



where 



3A 1 rrr, tk^t* — r 1 AN B(h,p) A / 1 B(h,p) 



(4.12) 



A(h,N) = ^ + R-2N~ 1 AN, B(h,p) = 2p 2 - 4p% + . (4.13) 

Gat o27r z G^ 

In terms of the Lagrangian variables, the constraint $4 = yields 



i I Gjv-t? 6A \ 6A / 

1 + 6A /V 

This is nothing but the Hamiltonian constraint in j2]. Similarly, plugging (14. lip into (I4.5p . 
the constraints $5 = yield 




1 + Gal 7? 

1 T 6A /v 



(4.15) 



These are the momentum constraints in [2]. We, however, note that it is more appropriate 
to regard the following linear combination 

'I'-,, = $ 5J - $idjN - $ 3 <9,e , (4.16) 

as the momentum constraints. This is because these are the constraints that generate the 
spatial diffeomorphisms for the phase-space variables, h^, p % \ N and e, rather than $ 5j - (see 
Appendix [AT) . 

By an explicit computation, one can show that $2j and $5^ are first-class constraints, 
forming the constraint algebra flA.52p - flA.57p . and the rest are second class (see Appendix 



|X]for details). Hence, the constraints are classified as 

#(phase-space variables) = j^[hij,p % \ N,p N , Ni,p 1 ^, e,p e ) = 2(6 + 1 + 3 + 1) = 22 
#(2nd-class constraints) = #($1, $3, $4, $ 6 ) = 1 + 1 + 1 + 1 = 4 , 
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#(lst-class constraints) = #(&2j, $5j) = 3 + 3 = 6. 
Consequently, the number of physical degrees of freedom (DOF) in n-DBI gravity reads 

DOF of graviton = - (22 - 4- 2x6) = 2 + 1. (4.17) 

Indeed, as advertised, we find an extra degree of freedom as compared to GR, which, as seen 
in the previous section, is a scalar mode, hence a scalar graviton. 

5 Potential Pathologies 

We shall now address some potential pathologies associated with the existence of the scalar 
graviton, which are known to afflict HL gravity. 

5.1 Vanishing lapse 

One of the most serious problems of HL gravity is the absence of dynamics discussed in 
[TO] . The problem is that the time flow of the Hamiltonian constraint yields an independent 
constraint on the lapse, which, as it turns out, requires the lapse to identically vanish in 
asymptotically flat space-times, implying that HL gravity does not have any dynamics. The 
authors of [10] emphasise that this problem is intimately related to the fact that in HL 
gravity the Hamiltonian constraint is second class. 

In n-DBI gravity, the Hamiltonian constraint is also second class. This raises the concern 
that this model might also lack dynamics. However, this problem is evaded in a similar way 
as the consistent extension of HL gravity does. Namely, the time flow of the Hamiltonian 
constraint, given by ( ]A.58j) in the full theory, does not yield an additional constraint unlike 
HL gravity; rather it determines the Lagrange multiplier A3, since the Hamiltonian constraint 
$4 depends on the auxiliary field e (and implicitly on the lapse N through e by solving 
$6 — 0) and thus {$3, $4} 7^ 0. This stems from the nonlinear lapse dependence in the n- 
DBI action and can be seen more clearly in the linearised theory in Section EXJjQ the time 
flow of the Hamiltonian constraint $4 is generated by its commutator with the Hamiltonian 
density (13.61) . Because of the lapse term A 2 n in the Hamiltonian constraint descended from 
the nonlinear lapse term (An) 2 ~ (iV -1 AiV) 2 in the action, the coefficient {$!,$ 4 } of the 
Lagrange multiplier Ai is non-vanishing. Thus, such time flow determines Ai rather than 
imposing an extra constraint on the lapse. 

To summarize, the nonlinear lapse dependence intrinsic to n-DBI gravity rescues the 
model from the absence of dynamics. 

7 In the full theory, the nonlinear lapse dependence is somewhat obscured by the introduction of the 
auxiliary field e which linearises the lapse dependence. 
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5.2 Short distance instability 



In [TT] it was found that the scalar graviton in HL gravity develops an exponential time 
growth at short distances in generic space-times in the linearised approximation. This sug- 
gests the presence of a universal short distance instability in HL gravity. Such problem, 
however, should be considered with some caution. The exponentially growing mode involves 
the lapse; but the lapse in asymptotically flat space-times is forced to vanish everywhere [10J, 
as mentioned in the previous section. If the lapse must collapse everywhere, there would be 
no way to develop an exponentially growing mode involving the lapse. Thus there appears 
to be a contradiction between these two issues. This might imply that the unstable scalar 
mode found in [11] fails to obey the boundary condition at spatial infinity, when extended to 
long distances. In any case, as we have just shown, the lapse does not vanish in our model 
and so it is sensible and important to study whether the scalar mode leads to any sort of 
instability or not. 

We shall now show that the analysis that reveals an exponential time growth at short 
scales in HL gravity, does not allow the same conclusion when applied to the scalar graviton 
in n-DBI. We begin by recasting the equations of motion of n-DBI gravity: 

(d t -C N )h lJ -2NK ij = , (5.1) 
R + K 2 - KijK ij = 2G N A c (e) + 2e~ 1 Ae , (5.2) 
V J (Kij - h tJ K) = -e-^e (Hy - h %J K) , (5.3) 
(Sit - C N ) (K ij - h ij K) + N (R ij - KK ij + 2K il K\ - h ij K mn K mn ) - V ViV + h ij AN 
= -e-\d t - C N )e (K ij - h ij K) + e~ l N (VV - h ij (V t lnN)V l ) e , (5.4) 

where e = (l + &l K)~* and the cosmological constant A c (e) = J£- - 1 - ^). The first 
equation defines the extrinsic curvature. The second and third equations are the Hamiltonian 
and momentum constraints, respectively. The last equation is the evolution equation [2]. 
In this form, the R.H.S. of ( 15.21) — ( T5~4|) represent the modification from GR. In fact, when 
the auxiliary field e is constant, the R.H.S. vanish except for the cosmological constant in 
the Hamiltonian constraint (15. 2p . and these equations become those of GR, as remarked in 
Section [2~T1 As we have seen in (I3.24p in Section l3.1.2[ the time-derivative of the Hamiltonian 
constraint succinctly captures the scalar mode dynamics, and it can be written as 

Ae - KNAe - (N&K + 2Kd l N) d { e - N^ANe = . (5.5) 

We now consider small perturbations around a generic background in the gauge = 0: 

= hij + 7y , (5.6) 
K ij = + K ij ; ( 5 7 ) 

N = N + n. (5.8) 

The details of the following analysis can be found in Appendix IB. 11 We assume that the 
background space-time varies over the characteristic scale L; R\j ~ ~ 1/ L, TZ ~ 
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1/L 2 , din N ~ 1/L, and e ~ 1 + 0(1/L 2 ). Since we are interested in a potential universal 
short distance instability, we focus on scales much shorter than L; 

w,p>l/L, (5.9) 

where the space-time is nearly flat. We can then neglect terms with space-time derivatives 
of the background fields relative to those of the perturbed fields. Assuming the Fourier 
decompositions 

llv K l \n^e- M ^ , (5.10) 

the perturbations of (JSHJ-O and ([53]) yield flR6|) - (lR9|) and ( lB14j) . respectively. We 
then find that for the scalar mode to be present 

u ~ Np ~ i (NK + ON) ~ - . (5.11) 

L 

At first sight, it might seem to imply the existence of an exponentially growing mode, 
suggesting an instability. This relation, however, violates the validity of our approximation 
( 15. 91) and thus cannot be satisfied. We conclude that this analysis does not reveal a potential 
universal short distance instability as discussed in the HL case [TT]. Note that this also 
implies that there is no scalar mode of the type (n, 7) ~ (n(cu,p), j(cj,p)) e lu ' t+ip ' x at short 
distances @ 

The same conclusion can be reached by an alternative analysis using the Stiickelberg field 
cf). We leave the details of the analysis in Appendix IB. 21 The equation of motion for <fi can 
be obtained from the action (I2.12p . The preferred choice of time-like vector n = (-N, 0, 0, 0) 
corresponds to <fi — t. Indeed, the equation of motion flB.20j) for <f> reduces, when </> = £, to 
the time derivative of the Hamiltonian constraint (I5.5p . 

We now consider perturbations in unitary gauge = t + ip and N l = 0, expanding (15. 5p 
to linear order in ip. To leading order, we find 

A 2 <p- (NK -d t \nN) A 2 <p + 2di\nN '&A<p = . (5.12) 

This indeed yields (15. lip found in the previous analysis. 

We conclude by remarking that, as advertised, <p{t,x) can be identified with the broken 
gauge parameter T(t,x). Since in flat space-time (I5.12p yields A 2 <p = 0, we obtain precisely 
the scalar graviton equation (I3.24p whose origin can be traced to the nonlinear lapse term 
(observe also that (I5.12p is the time derivative of the linearised Hamiltonian constraint). 
This provides a consistency check of our computations. 



8 As we have seen in section T3.1. 21 the scalar mode does exist, but it has ui = 0. This can also be seen in 
the Stiickelberg formalism, cf. (|5.12l) . The scalar mode is neither exponentially growing nor oscillating to 
linear order. We will discuss more about the scalar mode in Section [5] 
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5.3 Strong coupling 

Naively, n-DBI gravity admits the GR limit by tuning the parameter A — > oo and q — y 1 with 
X(q — 1) fixed finite. However, as is well-known, it is subtle and far from obvious whether 
the extra scalar mode actually decouples in the IR in the putative GR limit. A famous 
example is the vDVZ discontinuity of massive gravity [15] which can be attributed to the 
strongly coupled nature of a scalar mode, as elucidated in [IS] . Massive gravity has three 
extra modes; a transverse spin one and a scalar mode. As it turned out, the scalar mode 
remains strongly coupled in the massless limit and thus the GR limit does not exist. 

The original form of HL gravity has a similar strong coupling problem [Til D2] ; there is an 
energy scale A s above which the scalar mode self-coupling becomes strong. For HL gravity 
to flow to GR in the IR, the strong coupling scale A s needs to be sufficiently high so that 
the coupling becomes weak and virtually decouples in the IR. However, as it turned out, the 
naive estimate of the strong coupling scale yields a value too low, A s = \Ahl — 1 Mp, where 
Ahl is an anisotropy parameter and the coefficient of the K 2 term; this is most manifest in 
the Stiickelberg field formalism. The part of the Stiickelberg field action, coming from the 
K 2 term, is proportional to (Ahl ~ l)Mp. This is the most relevant part of the action in the 
IR. This implies that the coupling is g — 1/ (-\Ahl — 1 Mp) and thus irrelevant. Although 
the coupling becomes weak in the IR, the strong coupling scale A s is too low, approaching 
zero in the putative GR limit Ahl = 1. In fact, as argued in [TTJ, the situation appears to be 
even worse. The quadratic action in flat space lacks time derivatives, while there are time 
derivative interactions. This may imply that the high frequency modes are always strongly 
coupled. Indeed, one can refine the estimate of the strong coupling scale by making the space- 
time slightly curved (with the characteristic length scale L) which introduces a quadratic 
term with a time derivative. Then the refined strong coupling scale becomes anisotropic and 
is estimated to be A u] = L 

-l/4 A 3/4 

and A„ = L 

-3/4 A l/4 

[TTJ. In the flat space limit L — > oo, 

these are in fact zero. 

In n-DBI, the situation turns out to be more subtle. The Stiickelberg field action in flat 
space in unitary gauge yields 

- , (5.13) 

to cubic order, as shown in Appendix IB. 31 The quadratic action is invariant under the 
scaling E — > sE, t — > s~ 1 t, x — > s~ 1 x, and tp — > s _1 <£>. First note that the last interaction 
is irrelevant as compared to the middle ones. The higher order interactions have a similar 
structure. So the most "relevant" terms are of the middle type and actually (classically) 
marginal, with coupling A'V 2 , for the n-th order in ip% However, due to quantum effects, 
these operators may become either relevant or irrelevant. For the cubic terms with coupling 

9 Most easily seen by redefining 4> — > ^/\~<fi. 
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g 2 = g 2 Q l + 2c(? 2 ln 



M P 



(5.14) 



where A is the energy scale and the UV cutoff scale is Mp. Therefore, the strong coupling 
scale is 



This implies that in the putative GR limit A — > oo, the strong coupling scale is about 
the Planck scale and thus sufficiently high. However, for this coupling to be (marginally) 
irrelevant, the constant c must be negative. Given the fact that this is a scalar field theory 
and the higher derivative nature of the operators, it seems plausible that this is actually the 
case. Therefore, n-DBI gravity is likely to be free from the strong coupling problem. 

6 Discussions 

The first main result of this paper is to have established that n-DBI gravity has three degrees 
of freedom: the usual two tensorial modes of the graviton plus a scalar mode which, however, 
is not a propagating particle with a definite dispersion relation. The second result is that the 
arguments for three pathological features of the original formulation of HL gravity do not 
carry over to n-DBI gravity; these arguments concern the absence of dynamics, a universal 
short distance instability, and a strong coupling issue. The reason why n-DBI looks healthier 
than HL gravity, in these respects, has a common ground: the nonlinear dependence on 
(spatial derivatives of) the lapse in the action. 

The fact that the same arguments for the pathologies in HL gravity do not directly apply 
to n-DBI gravity does not mean, of course, that the latter is free of pathologies. Although 
the arguments presented herein are reassuring, in order to establish an healthy behaviour, 
n-DBI gravity must be further studied. In particular, we found that the scalar mode grows 
linearly in time in the linear approximation and is thus at the threshold between stability 
(oscillation) and instability (exponential growth). In other words, the scalar mode can be 
either marginally stable or unstable!^] The (in) stability depends on the details of the scalar 
mode interactions. We thus need to extend the analysis performed in section 13.1.21 beyond 
linear order to establish the (in) stability of the scalar mode. 

To study the effect of nonlinear interactions, we find it useful to work in the Einstein 
frame (12 .7p . Note that the background space-time remains flat under the frame change. A 

10 Here we assume that the IR divergences can be properly regularized. 

11 As an illustration, consider a simple mechanical model. For canonical kinetic term, a linear time growth 
is the behaviour seen in a flat potential. The flat potential does not indicate an instability and is marginally 
stable. However, any potentials which are flat to quadratic order all lead to the linear time growth for per- 
turbations in the linearised approximation. The simplest examples are cubic and (convex) quartic potentials. 
The former is clearly unstable (inverse-squared blow-up), whereas the latter is stable (oscillation). 




(5.15) 
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Figure 1: The potential V(x) defined in (12.81) with q = 1 

key observation is that, as suggested in Appendix (CJ the auxiliary field x can essentially 
be thought of as the time derivative of the scalar mode, i.e., x ~ T ~ if. In other words, 
the proof of the stability amounts to showing the stability of the x fluctuation. In fact, the 
equation of motion for x i n the N l = gauge 

d t (N~ 2 X ) + A X - ^(N-'K) - h/'{x) = (6.1) 

seems to indicate the stability against small perturbations around x = 0; as is clear, the 
potential V(x) plotted in Fig{T]is very stable around x — 0- Note, however, that the lapse 
N and the spatial metric also furnish the scalar mode. Moreover, the auxiliary field x is 
further constrained by the other equations of motion. So, strictly speaking, we need to deal 
with the coupled system of x, N and hij (in the N l = gauge). There is, however, a shortcut: 
it would be sufficient to show that the solutions of (16. ip for x i n generic backgrounds are 
oscillatory. 

The force in (16. ip consists of the strongly attractive force ^N 2 V'(x), the repulsive force 
— N 2 Ax = N 2 p 2 x, the friction d t \nN 2 x ~ L^x, and the external force |JV 2 9 t (iV -:l iif) ~ 
L~ 2 . The last three forces may work against stability. These become more dominant when x 
is smaller. However, the small x behaviour is well approximated by the linear perturbation. 
As shown in Appendix [Cj x is marginally stable to linear order. Meanwhile, for larger x, the 
attractive force, being exponential, may quickly become dominant, as x increases, and pull 
it back to smaller x- Hence it seems plausible that \ oscillates about x = 0- We therefore 
conclude that n-DBI gravity is likely to be stable against the scalar mode perturbations. 

To close, we have considerably improved our understanding of the scalar mode and pre- 
sented arguments for its healthy behaviour. However, the scalar mode still remains somewhat 
elusive; it does not appear to behave as a conventional particle with a (non-)relativistic dis- 
persion relation around near-flat space-times. So it is rather unclear what signatures exactly 
it might leave as an observable. On this score, however, we anticipate that the scalar mode 
could be the source of scalar perturbations in the Cosmic Microwave Background Radiation, 
when applied to inflationary cosmology pQ, and we hope to return to this question in the 
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near future. Finally, we believe that models such as n-DBI gravity are interesting learn- 
ing grounds for deepening our understanding of (non-)relativistic gravity. Needless to say, 
hurdles are high to self-consistent ly modify Einstein's elegant theory of general relativity. 
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A Computations of constraint algebra 

In this Appendix, we give an explicit computation of the constraint algebra and classification 
of the constraints presented in Section HI To facilitate the computation, we introduce smooth 
test vector fields, £ M = and rf = (rf,rf), which fall off fast enough to suppress all 

the boundary contributions [T7]. Henceforth, we define the smeared constraints: 




($ l 2 V,iW +V j ($ 2i W))(x) , (A.9) 



(A.10) 



(A.8) 



(A.7) 



(A.5) 



(A.6) 



(A.4) 



(A.3) 



(A.2) 



(A.l) 
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where $1, $2«, $3, $4, ®5i, ^6 and $5; are the constraints defined in Section HI The idea is 
to compute the commutators like 



(A.ll) 



{$ 4 (£ U ),<W)} = J rfVX u (y)^(x){$ 4 (y),$ 5j (x)} 

and read off the algebra of local constraints from the R.H.S. 
The basic non-vanishing Poisson's brackets are given by 

{p l Ky),h k i{x)}= l -{5l5l+5\5i)5{y-x) , 

{p N (y),N(x)} =8(y-x) , 
{p%{v),N j (x)} = S t j S{y-x) , 
{Pe(y),e(x)} = S(y-x) . 



(A.12) 

(A.13) 
(A.14) 
(A.15) 



To compute the Poisson brakets of constraints, we take the variations of the smeared con- 
straints: 



5p N 
s Pm 

Sp e 



(A.16) 

(A.17) 
(A.18) 



K 



Z h mn^ 4 _ __ h mn _ _ p 



+ -VhR mn - ^(V m V"e) 

K 



2n f 

Wl v 



p ml p n _ _ p mn p 



Sp mn 
be 

bh mn 

SMC) 
5p mn 

jge(g) 
5h 



+ 



e 

K 



K 



2k ( 



-(V m V n e°)e + h mn {V a e){Va£ ) + h mn e{A£°) 
1 



K 



Pmn 2 h m nP 



P ab pab--P')+—R + 



e 2 Vh 



Vh„ QXVh 



K 



KG 



N 



1 

e- 



K 



(A.19) 
(A.20) 
(A.21) 
(A.22) 
(A.23) 



~-h mn <S> 6 + 



Nk 



h mn p ab Pab --P 2 )- 



2kN ( 
e 2 Vh V 



P Pi - 2? P 



K K 
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+ 



Vh 



-N(V m V n - h mn A)£° + h mn {V l N){V£°) 



= -e 
= e 



5p mn 

5N 

jgg(g) 
be 

^S? } (e) 

5N 
5p N 
5e 



2Nk f 1 

^"y= \Pmn - ^h mn p 



K 



P ab p ab -=:P Z ) +—R + 



Vh„ Q\Vhf_ 1 



( ab 1 2 



12\nVK 



\ \ 1 e 2 



5N> 

jv 

^Pn 



= di(Cp N ) , 

= -e^iV = -£ ( N , 

= ~CdiC =-£{:€, 



Using these, it is tedious but straightforward to compute the constraint algebra: 



{Si(l/),* 6 (aO} = 



) + —R + 



K 



kG 



N 



5(y - x) 



+ 2— A% - x) , 

K 



{^(y), $ ( $(x)} = Q 1 (x)d yi 6(y-x) , 



{cJ> 3 ( y ),<D 4 (x)} = - 



P° V - ) + — # + 



ft 



ftG 



A? 



1 - 



+ A% - x) , 

ft 

2Nft / h 1 9 \ 12XNVK' 
P Pab - ~P 



{*-<*).*■<*)>= 2 ,, ^ 

{*3(v),S^(z)} = M*)M(!/-z) , 

{$ 4 (y), $4(a:)} = $5(2/)^% - x) - $^(x)^5(y - x) 

- ^d y ked yk 5{y - x) + ^d xk ed Xk 5(y - x) , 



5(y-x) , 



(A.24) 
(A.25) 
(A.26) 
(A.27) 

(A.28) 
(A.29) 
(A.30) 
(A.31) 

(A.32) 
(A.33) 



(A.34) 
(A.35) 



(A.36) 

(A.37) 
(A.38) 

(A.39) 



21 



{$ 4 (y), $5j(x)} = <S>i{x)d yj 8(y - x 
2N Pl 



^,{x)5{y-x)-2^d xi e[A 



AN\ 



N 



J S(y-x) , 
(A.40) 



{$ 4 (y),$ 6 (x)} 
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V m V m e V m V n iV 
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d Xm S(y-x) 
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S(y - x) 



{$ 5j (y), $5i(x)} = <f> 5j (x)d y id(y - x) - <$> 5i (y)d xj 6(y - x) , 
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V e 3 



5(y - x) 



(A.41) 

(A.42) 
(A.43) 
8{y - x) 

(A.44) 



{My),M*)} = -dv 
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V e 3 
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{*g°(y), CV)} = ^(i)^*(y - x) - *S"(!/)3tf*(v - s) 



AG) 



(A.45) 
(A.46) 

(A.47) 

(A.48) 
(A.49) 



{$!(y),$ 2j (x)} = {$ 2i (y),$ 2j (x)} = {$ 3 (y),$ 2j (x)} = {$ 4 (y), $ 2j (x)} = {$ 5l (y), $ 2j (x)} 
= {$ 6 (y),$ 2j (x)} = {<&$(y),<M*)} = {$?(y),$ 2j (x)} = . (A.50) 



In order to classify the class of the constraints, it is more appropriate to choose 

($l,$ 2i ,$ 3 ,$4,$5j,$6) 



(A.51) 



as a set of independent constraints. $5^ is defined in flA.101) and given by a linear combination 
$ 5 i — diN&i — <9je$ 3 . As is clear from the above computation, <3> 2j - and $ 5 j commute with 
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all the constraints: 



{*i(y), 


= {$ 


2i (j/),$ 2i (x)} 


= {$3(j/),^(x)} = {$ 4 (y),^(x)} 






= {$ 


5i(y),$2j(x)} 


= {$ 6 ( 2/ ),$ 2 ,(a;)} = 0. 


(A.52) 
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{$ 3 (y),$ 5j (x)} 
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(A.55) 


{$ 5i (y),$ 5j (x)} 
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) - $ 5j {y)d x id(y - x) pa , 


(A.56) 


{* 6 (v), 


= $6< 


x)d y jS(y - x) 


pa . 


(A.57) 



It is easy to show that the set (1A.51I) is complete. Namely, the time flows of the secondary 
constraints do not give rise to any new constraints: 

$ 4 (x)= / d 3 y{H e j$ il (y),Mx)}+ E / *v{*a(v),H*)}\>+ [ d*y{%{y)^{x)}\ 2i 

J o=l,3 ^ ^ 

pa , (A.58) 
*sj(x)= [ d 3 y{n e JX(yl^)}+J2 I d?y{^ a {y)^{x)}\ a + [ d 3 y{¥ 2 (y),$ 5j {x)}\ 2 ,, 



a=l,3" 

pa , (A.59) 

$ 6 (X)= / A{K ( DBl(Z/)^6(x)}+ £ / A{$«(y)^6(x)}A a + / A{$ 2 (l/), $ 6 (x)}A 2i 

a=l,3^ ^ 

pa . (A.60) 
Since the Hamiltonian density takes the form 

U^S m = - + N j $ Bj ) - -Vh(eAN - NAe) , (A.61) 

one can see that (1A.58I) and (1A.60I) determine Ai and A3, while (1A.59|) is trivially satisfied. 
Hence, there are no additional constraints from these time flows. Having established the 
completeness of the set (1A.51[) . we conclude that $ 2 j and $ 5 j are first class, and the rest are 
second class. 

Finally, let us end this appendix with a comment on the generator of the spatial diffeo- 
morphism. The generator of the spatial diffeomorphism acts on a phase-space variable 
A as 

{A(y),g(?)} = £eMv) • (A-62) 
Since Pn,P 1 ^, and p e are primary constraints, we only need to consider the reduced set 
of phase-space variables, (hij,p l \ N, N,e). The spatial diffeomorphisms for this set are 
generated by ($ 5 , — $ 5 , ^e^)- They are indeed all generated by the first class 

constraints, as can be seen from flA.7p - flA.10p . 
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B Computational details of perturbations 



In this appendix, we show some details of the computations in the perturbative analysis of 
the scalar mode in Sections 15.21 and 161 



B.l Perturbation of the equations of motion 

The linearised version of equations ( 15.ip -( |5~4l) and (15. 51) can be obtained in a fashion similar 
to [UJ. To the approximation explained in Section [5721 we find, in the gauge N l = 0, 

7y - 2NK tj - 2K ijn = , (B.l) 

C 

2KijK ij - 2kK - V^^ij + A 7 = -^-Aa , (B.2) 

oA 

yj K .. _ ViK _ Z Rjk v . ljk + K^V^ij + ~#iiV J 7 = ^ [K l3 - htjK) V ] a , (B.3) 

k ij + f j K - h ij k - j k ih ij K kl - V + h ij An 

N C 

+ 17 (V fc VV* + V fc VY fc - - W 7 ) = -T^iVVVa , (B.4) 
2 12 A 



and 

Ad = Ad t 



W'to- - A 7 + 2KijK ij - 2kK - 2N- 1 ^An - Va ij VjN + -V^ViiV 



= . 
(B.5) 



The L.H.S. of f]B.l|) ( ]BT4|) are the same as those in GR and agree with the A = £ = 1 case 
of [11]. Eq. flB.5j) . however, is very different. 

We are interested in perturbations of the modes with wavelengths much shorter than the 
characteristic scale L of the background, i.e., u,p ^> 1/L, where the space-time is virtually 
flat. First, (IB.5|) enforces a to the constant. Thus the R.H.S. of (IB.ll) - flB.4l) are always 
negligible and thus we have in Fourier space 

iuY j + 2(Nn ij + K ij n) = , (B.6) 
( PiPj - p 2 6 lJ ) 1 lJ - 2(Kk - K ijK ij ) = , (B.7) 

p . K V - p i K - Ik^p* + K jh p k ^ + - 1 K ,J ,>f = , (B.8) 
iu(5 ij K - K ij ) - {p 2 5 ij - pY)n - iu{Ry j - 5 ij K k a H ) 

- j(PkP J l tk + p k py k - P 2 1 IJ - pV't) = • (B.9) 

Plugging (IB. 61) into ( 1B.9I) and discarding sub- leading terms linear in ~ O (1/L), we 
simply obtain the fluctuation equation of GR in flat space-time: 

^(S ij K-K ij )+iC0 ^ - n + N 2 (^K ij + ^ (pV'k - p kP >K ik - Pkp'K^ = . (B.10) 
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Contraction with 5y and pj, respectively, yields 

iuK-p 2 n=0, (B.ll) 
p^K — pjK^ —0 . (B.12) 

Using these relations, one can find that ( IB. 101) gives the massless dispersion relation 



uj 1 = NY . (B.13) 

In this approximation, the rest of the equations, the Hamiltonian and momentum constraints 
(IB.7j) and (1B.8|) . are automatically satisfied. In GR, we can set n = by using the residual 
gauge symmetry. Thus we are left with a massless transverse traceless tensor mode, i.e., the 
usual graviton. 

In n-DBI gravity, we do not have liberty to gauge away n. However, we have ( IB. 51) to 
take into account. Using (1B.6|) . it becomes 



[iuN (Kij - K5 ij ) + N 2 (p 2 6 ij - p iPj + d.^h), In N - 2ip i d j In N)] 

+ n [-iup 2 + NKy (p 2 5 ij - p { pj + ip k d k In NSy - 2ip i d j hxN)] = . (B.14) 

To leading order, this yields 

N 2 K tj (p 2 5 ij - pip^ - inujp 2 = OP nup 2 = . (B.15) 

Since our approximation is valid only for u,p 3> 1/L, this implies that n = and we are 
again left with a graviton. To the next order, however, this reads 

n ij [iuN (Kij - K5ij) + N 2 (5ijip k d k miV - 2ip i 6> j InJV")] + n[-iup 2 + - PiPj)] = 

(B.16) 



Due to (IB. Ill) , the k u terms contribute. For n to be non- vanishing, we must have the relation 



u~Np~i(NK + NdhiN) ~ 1 . (B.17) 

However, once again, this cannot be satisfied. Hence there is no scalar mode of the type 
(71,7) ~ (n(co,p),j(cj,p)) e lU}t+w ' x with \u>\, \p\ 3> l/L, including the one with imaginary uj. 

B.2 Perturbation of the Stiickelberg field 

To obtain the equation of motion for the Stiickelberg field 0, we vary the action (12. 12ft or 

(EH 

5S ~ e5IC ~ e5 (D a (n a D n fi )) ~ (D Q (n^e) - ifD a e) 5n a , (B.18) 
where e = (l + ((% + /C))~* and /C = -2D a (n a D p n p ). Using 



*n« = ^^ + ^ — , (B.19) 
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we find 

nQ ( ( 9 <x + n a n lI d )(n p dpe) - D$nP{d a + n a n a d a )e\ 

D ^ — j = . (B.20) 

It can be shown that this becomes (15.51) when n = (-N, 0,0,0), i.e., = t. Note that 
this equation is invariant under <fi — > /(</>), as it should be, owing to the fact that n a (d a + 
n a n a d a ) = 0. This also implies that the quantity in the parenthesis is proportional to the 
space-like vector n^Dpria tangential to the hypersurface. 

Now we consider perturbations <fi — (j) + (p, expanding (IB.20j) to linear order in (p. We 
work in unitary gauge, <f) = t and N l = and use 

/ , \ -1/2 

Sn a = (0,N&(p) , 8 [V^X) = -N<p , 5K = N- 1 V l (N 2 d t p) , (B.21) 
as well as 

Se = -^re 3 5!C , -~<WC = (2K - ^^(JV^V^JV 2 ^^)) + N&ydiK . (B.22) 
12 A 2 



The terms coming from perturbing n a and 1 / y—X in (1B.20j) contain at most 3 ^-derivatives 



while those from perturbing e contain 6 and 5 ^-derivatives. Thus we only need to consider 
the latter. It also suffices to keep 3 and 2 5e-derivatives. Then we find 

A5e + KNA5e = . (B.23) 

In terms of ip the leading terms (with 6 and 5 derivatives) yield (I5.12p . 

B.3 Perturbation of the Stiickelberg field - 2nd order 

To find the cubic action of the fluctuation p in flat space-time (with q — 1), we expand the 
action as 



S = 1 / ,] [ r 
IQttGn 



4 V 6A / 8 V 6A 



(B.24) 



The first term is a surface term and does not contribute to the equation of motion. In unitary 
gauge, the time-like vector takes the form 

n = , rii = . (B.25) 

y/il + ff-ftip)* ^(l + vi) 2 -^) 2 

This can be expanded as 

- n = 1 + i(«9^) 2 - + 0(<p 4 ) , (B.26) 

- n t = dup - pd t p + p 2 d^ + ^pf + 0{p> 4 ) . (B.27) 
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In the flat background we have 



K, = -2 [fib (n°d n°) + di (n i d j n j ) + d (n ^) + d { (n^n )] , (B.28) 

and we find to quadratic order 

/C = 2[A^ + fi (^A^)-fi i (S^A^)] . (B.29) 

Hence the cubic scalar field action is given by 
1 



192ttA 



d x 



1 1 

-(A0) 2 - (0Ay3 - fi^Afi^) + ^(A0) 3 



(B.30) 



B.4 A nonlinear analysis of the Stiickelberg field 



The scalar mode ip obeys the equation of motion (1B.20j) . To determine whether the scalar 
mode leads to an instability or not, we need to study (IB. 201) beyond linear order approxima- 
tion. As we discussed, to linear order, the equation of motion is simply 



A 2 if = 



(B.31) 



and the solution is 

(p(t,x) = e [<f (x) + (fi{x)t] , (B.32) 

where ipo(x) and ipi(x) are arbitrary functions of space. We have included the factor of e <C 1 
for later convenience. The fully nonlinear solution can in principle be found systematically 
order by order in e expansions: 



cp{t,x) = e [<p Q (x) + <fi{x)t) + ^ e n (p n (t, 



x) 



(B.33) 



n=2 



The higher order fluctuations ip n (t, x)'s are determined in terms of the initial data (ip (x), <fi(x)) 
and of order n in powers of (spatial derivatives of) fo(x) and f\{x) and polynomial in time 
t. Using (IB. 301) . for example, the next-to- leading order fluctuation can be found as 



ip 2 (t,x) = ^\x)t 2 + ^tp¥\x)t 3 



(B.34) 



where 



AVf =A [2A^oA^! + 2d i ip Ad i ip 1 + Afi^o^i] - (A^AVi + Afi^oAfi^i) , 
A*<p® =A [2(A^) 2 + 2(fi^ 1 ) 2 + dMAdnpi] - (AViA^x + d m Ad lVl ) . (B.35) 

Clearly, the late time behaviour of the scalar mode requires the knowledge of all order 
fluctuations. Thus, to see whether the scalar mode yields an instability or not, we need to 
resum the infinite series (1B.33j) . However, this seems to be out of our reach and we will 
instead resort to an alternative analysis working in the Einstein frame. 
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C Linear perturbations in the Einstein frame 

We consider linear perturbations of the scalar fields around flat space-time in the Einstein 
frame (12.71) . The scalar field Lagrangian density reads 

AkGnC*^ =2^ 2 + 4iP(E + AB) - (An - 2ip)Aip 

+ A X (e + AB + 2V>) + 6x 2 - 6 X A X - ^X 2 • (C.l) 

This is the Einstein frame counterpart of (13.31) . The equations of motion are given by 

4> = -X , (C2) 

Ai, = -A X , (C.3) 

E + AB + An = -x , (C.4) 

Aip = , (C.5) 

which clearly reduce to those of GR for constant Xi P ms 

1 4A 

X + A X + -{E + AB + 2ijj) + —x = . (C.6) 

This is the linearisation of (16. ip . In the E = gauge the general solution can be found as 

B =B (x) + B 1 (x)t , (C.7) 

n = -B 1 (x), (C.8) 

i> =0 , (C.9) 

X=Xo{x), (CIO) 



with Xo( x ) related to Bi(x) by 



1 9 \ 

AB 1 = -3A X o- T ^Xo ■ (Cli; 



G 



N 



Here we have again imposed the boundary condition that all the fields fall off at spatial 
infinity. Note that Xoi x ) is essentially B\(x) which is the degree of freedom responsible 
for the linear time growth. This suggests the identification x{t, x ) ~ T(t,x), that is, the 
auxiliary field x can be regarded as the time derivative of the scalar mode. 
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